<N 



CO 

O 



LPSC 1084 

DCPT/10/90 

IPPP/10/45 



On SUSY GUTs with a degenerate 
Higgs mass matrix 

F. Briimmer^, S. Fichet^, S. Kraml^, R. K. Singh ^ 

^ Institute for Particle Physics Phenomenology^ Durham University^ 

Durham DHl 3LE, UK 

O ^ Laboratoire de Physique Subatomique et de Cosmologie, UJF Grenoble 1^ 



CNRS/IN2P3, 



^ 53 Avenue des Martyrs, F-38026 Grenoble, France 

Q^ ^ Institut filr Theoretische Physik und Astronomic, Universitdt Wilrzburg, 

D-97074 Wiirzburg, Germany 

Oh 

Oh 

(D 

l^^ Abstract 

^^ Certain supersymmetric grand unified models predict that the coefficients 

^_^ of the quadratic terms in the MSSM Higgs potential, 772^2 — ^Hi2 ^ 1^1^ 

CN and m^ = B^, should be degenerate at the grand-unified scale. We discuss 

some examples for such models, and we analyse the implications of this 

peculiar condition of a GUT-scale degenerate Higgs mass matrix for low- 

(^ scale MSSM phenomenology. To this end we explore the parameter space 

^^ which is consistent with existing experimental constraints by means of a 

Markov Chain Monte Carlo analysis. 
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c3 1 Introduction 

The quadratic part of the Higgs potential in the Minimal Supersymmetric Standard 
Model (MSSM) may be written as 

V = ml\Hi\^ + ml\H2\^ + ml{H2Hi + h.c). (1) 

Here ml 2 are given in terms of the soft SUSY breaking masses 7tz|^^ ^ and the super- 
symmetric Higgsino mass /x as m^ 2 = ^Hi 2 + I/^P- The soft parameter ml is often 
also called 5^; the phases are defined such that mg > at the electroweak scale. 
Hi 2 are the lowest components of the down-type and up- type Higgs superfields 
(which we will also denote by Hi^2)- 



In models of gravity-mediated SUSY breaking, the Higgs mass parameters are 
usually generated at the GUT scale or the Planck scale by some mechanism which 
breaks supersymmetry. (Even the /x parameter, although it preserves supersymme- 
try, should be generated by the SUSY-breaking mechanism in order to explain why 
it is of the same order of magnitude as the other Higgs mass parameters.) They 
should then be evolved down to the electroweak scale according to their renormal- 
ization group equations. At the scale where the Higgs potential is minimized, they 
should satisfy the well-known inequalities 

ml ml- ml < 0, mj + mj- 2ml > 0- (2) 

The first inequality ensures that electroweak symmetry is broken, and the second 
one guarantees that the Higgs potential is bounded from below even in those 
directions in field space where the quart ic potential vanishes. 

An interesting property of certain UV-scale models is the relation 

m^l = ml = img, (3) 

holding at the UV scale, where the SUSY-breaking terms are generated. This 
is the defining relation for the class of models we are interested in, models with 
a degenerate Higgs mass matrix (DHMM). It is typically encountered in models 
where both MSSM Higgs doublets originate from a single chiral adjoint $ of the 
GUT group G, by a decomposition into Standard Model representations according 
to 

Ad(G) ^(l,2)_i/2 0(1, 2)1/2 0... .^. 

^ ^ Hi //2 . . . 

Suppose that there is some additional structure ensuring that, while the real com- 
ponents $''' + $ can acquire a tree-level mass from SUSY breaking, the Higgs field 
components of the imaginary parts $t _ $ remain massless. For instance, they 
could be the (pseudo)-Goldstone bosons of a spontaneously broken (approximate) 
global symmetry fllj2| or a mass term could be forbidden by higher-dimensional 
gauge invariance [3[|4]. Then the quadratic Higgs potential will be of the form 

V = m\Hi + H2){Hi + :^). (5) 

Thus the relation Eq. ^ obviously holds. A subtlety lies in the definition of the 
sign of TTZg, which may need to be changed by a field redefinition Hi ^ —Hi to 
ensure that mg > after running to the electroweak scale. 

In realistic models in which the UV-scale equalities Eq. ^ apply, renormal- 
ization group running should turn them into the IR-scale inequalities Eq. ([2]). At 
the same time the correct Standard Model couplings and masses should be repro- 
duced, and the spectrum of sparticle and Higgs masses should not be in confiict 
with experimental bounds. Whether this is possible for a given UV-scale model 
can only be found out by a numerical analysis. As discussed in |3||4], it crucially 
depends on the parameters of the sfermion sector. 

In a previous paper |4| a particular subclass of such models was investigated 
in detail with respect to their phenomenological prospects. More precisely, it was 



shown that in five-dimensional orbifold GUTs with gauge-Higgs unification and 
radion-mediated SUSY breaking, fuUy reahstic MSSM spectra can be found. Such 
models are well motivated as anisotropic limits of heterotic string compactifica- 
tions. 

Since the class of DHMM models is in fact quite large and diverse, it is clearly 
of interest to study in detail under what conditions on the soft terms realistic 
MSSM vacua and a viable low-scale phenomenology can result, and what are the 
consequences for experiment. This is the purpose of this paper. In Section 2 
we first discuss classes of SUSY GUTs which imply DHMM boundary conditions. 
In Section 3 we elaborate on RGE running and resulting soft term patterns. In 
Section 4 we then present results of a Markov Chain Monte Carlo (MCMC) analysis 
of two variants of DHMM models. Finally Section 5 contains our conclusions. 



2 Models 



2.1 Gauge-Higgs unification in 5d 

An interesting class of example models in which Eq. ^ holds is given by 5d orbifold 
GUTs with gauge-Higgs unification |3]-[6]. The fifth dimension is compactified on 
an interval whose size is given by the inverse GUT scale, and the GUT group 
is broken to the MSSM gauge group by boundary conditions. In terms of 4d 
superfields, the 5d gauge multiplet decomposes into a 4d gauge superfield V = 
—A^a^OO + . . . (/x = 0, . . . , 3) and a chiral adjoint $ = S + iA^ + . . . (where we 
have only written the leading terms in the ^-expansion, and V is in Wess-Zumino 
gauge). $ contains the MSSM Higgs fields as in Eq. Q. We can now choose a 
Kahler-Weyl frame such that the superpotential is independent of $ when setting 
the MSSM matter fields to zero. By 5d gauge invariance, the Kahler potential 
can then only depend on the combination $''' + $ on the quadratic level. The 
orthogonal combination $''' — $ ^ A5, being a 5d gauge field, is protected from 
getting a mass term. 

This can be seen explicitly as follows: Suppose for the moment that the gauge 
symmetry were just U(l). The action is invariant under 5d gauge transformations 

V ^ V + K + K, $ ^ $ + a5A. (6) 

Here A is an x^-dependent chiral superfield. The inhomogeneous transformation 
behaviour of $ shows that $ cannot appear in the superpotential if W is to be 
5d gauge-invariant, when setting the MSSM matter fields to zero. That is to 
say, it is always possible to shift harmonic terms from the Kahler potential into 
the superpotential, and any terms from W into X, but a particularly natural 
formulation is one where W and K are separately 5d gauge invariant. Consequently 
$ cannot appear in W (except in combination with other charged fields such 
as matter fields, which only give rise to Yukawa terms irrelevant to the Higgs 
potential, or light exotics, which we assume to be absent). 

The crucial observation is now [4J that in this manifestly 5d gauge-invariant 



formulation, the ^-dependent part of K must be a function of the unique gauge- 
invariant combination 

$ + $ - d^V. (7) 

This combination reduces to $ + $ on the zero-mode level. In other words, if there 
is no linear term in $, the low-energy effective Kahler potential for the zero modes 
has the structure 

x = /c(z^z') +:P(z^z') ($ + $)' + ... (s) 

Here the Z^ denote collectively the compactification moduli and general hidden 
sector fields. K cannot depend on the orthogonal combination ^45 = Im$ = ($ — 
$)/2 essentially because the transformation law for Im$ involves ImA, whereas 
the transformation law for V only involves Re A, and therefore the gauge variation 
of Im$ cannot be cancelled. 

A similar structure is encountered in realistic models. The gauge symmetry 
should of course contain the Standard Model gauge group, and $ should contain 
the MSSM Higgs superfields, so the abelian example is too simple. In a non-abelian 
model, the 5d gauge transformations read 

e^ ^ e^'e^e^, $ ^ e-^{d^ + ^)e^. (9) 

Gauge-invariant operators involving $ can be constructed from the covariant 
derivative V5 = ^5 + $ [T] . In particular the operator 

- e-^Vse^ = $ + ^^^ - d^V + (commutators) (10) 

(where V5 acts on e^ as V^e^ = d^e^ — ^^e^ — e^^) is the appropriate non-abelian 
generalization of ^ . Note that it is not gauge invariant by itself but transforms 
analogously to a field strength superfield: 

e-^Vse^ ^ e-^e-^Vse^'e^. (11) 

The lowest-order gauge invariant operator one can construct is in fact |7J 

tr (e-^Vse^)^ = tr ($ + ^^^ + (terms involving V) (12) 



since tr (e V^e ) vanishes identically, as can be seen from Eq. (10). 



As in the abelian case, any ^/-independent terms cannot depend on the or- 
thogonal combination $ — $t since it transforms as 

$ - $t ^ e-He^ - h.c. + ^5 (A - A^) , (13) 

while the gauge field transforms as 

V -^ y + A + A''' + (terms involving V). (14) 

Therefore there is no function of V whose gauge variation can cancel the inhomo- 



geneous piece in Eq. (13). 



We conclude that again W is ^-independent, and that K has the structure 

K = JC [Z\'Z'^ + y [Z\'Z'^ tr ($ + $t)' + . . . (15) 

The resulting quadratic Lagrangian for the zero modes of $ can then be written 
as 

C^neA= [d^Oipipy{z\Z'^ (Hi + H2){Hr+Th). (16) 

Here ip = 1 + F^'^d'^ is the conformal compensator of 4d supergravity, a non- 
dynamical chiral superfield whose F-term is the scalar auxiliary field of the 4d 
gravitational multiplet; its expectation value parameterises SUSY breaking in the 
4d gravitational background. A non-vanishing F'^ or non-vanishing F* will give 
rise to an effective Higgs mass matrix satisfying the relations Eq. (JSJ), with mass 
parameterq^ 



,2 



"^ifi = ^H2 



F'F'—^-\ogy, 
dZ'dZ' 



±li = F^ + F'-^logy, (17) 

dZ 



±B^ 



F^' + F'^My 



—- d^ 



The relations Eq. ^ are a direct consequence of 5d gauge symmetry, which 
is however not a symmetry of the 4d effective theory (but instead mixes the KK 
modes). Eq. ^ is therefore valid at the compactification scale, but it will be 
modified by radiative corrections below this scale. This is indeed necessary in order 
for the conditions Eq. ^ to be satisfied as strict inequalities, after renormalization 
group running to the electroweak scale. 

Models of this kind can be augmented by bulk hypermultiplets whose zero 
modes give the MSSM matter fields, and with brane fields with appropriate su- 
perpotentials to decouple unwanted exotics |5|. 4d Yukawa couplings are obtained 
from 5d gauge couplings, with their precise values controlled by the localization 
properties of the zero-mode wave functions, which in turn are tunable through 5d 
mass terms. The full model has the massless spectrum of the MSSM. With the 
additional assumption that SUSY breaking mediation is dominated by the F-term 
of the radion modulus, it can give rise to realistic phenomenology [4|. See Fig. [l] 
for a sketch of this kind of model. 

By similar arguments as above, the relations Eq. ^ also apply in a large 
class of heterotic string orbifold models [8]-[ll| with gauge-Higgs unification, if 
their moduli space admits a corresponding 5d orbifold GUT limit (for details see 



^As already alluded to, the signs of /i and m^ can be simultaneously flipped by a superfield 
redefinition Hi -^ —Hi (together with a corresponding redefinition of the down-type matter 
superfields to keep the Yukawa couplings intact). We adopt the usual convention that m^ > at 
the weak scale. Since the sign of m'^ can change during its renormalization group evolution, this 
implies that the proper high-scale sign in Eq. ( pT] ) is only fixed after specifying the rest of the 
model and tracking its RG running. The relation Eq. ([3| thus can be written as ml = ml = en ^3 
with ej:/ = ±1 to be determined accordingly. 
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x^ = brane x^ = ttR brane 

Figure 1: An unwarped model with gauge-Higgs unification. The field $ containing 
the Higgs fields and the 4d gauge fields (green dashed curve) are bulk fields with a 
flat profile. Third generation matter fields (red dot-dashed curve) are slightly localized 
towards the x^ = ttR brane. Matter fields of the first two generations (blue dotted curve) 
are effectively confined to the brane. 



e.g. I4lh2l). This is regardless of whether or not this anisotripic limit is actually 
realized at the point where the moduli are stabilized. The structure enforced by 
higher-dimensional gauge invariance persists independently. 

Radion mediation (corresponding to "modulus domination" in the above- 
mentioned string models) is a simple and elegant possibility to parameterise SUSY 
breaking in such models. However, in general there may be other contributions, 
in particular from brane-localized fields (see e.g. [3]). This allows for more gen- 
eral patterns of soft masses than the ones considered in |4], providing a strong 
motivation for the more general parameter space scan which we perform in this 
paper. 

2.2 Holographic GUTs 

A somewhat different example is the holographic GUT model of Nomura, Poland 
and Tweedie [13j. This model may be described in the "gravity picture" as a 5d 
theory on a slice of AdSs space between two 4d branes, a Planck brane (or UV 
brane) and a GUT brane (or IR brane). The bulk gauge symmetry is SU(6). The 
MSSM Higgs fields arise from a GUT-brane chiral superfield $ in the adjoint jj 
SU(6) is spontaneously broken to SU(4) x SU(2) x U(l) by the $ superpotential 
on the GUT brane, and explicitly broken to SU(5) x U(l) by boundary conditions 
on the Planck brane. In the 4d effective theory the gauge group is then given by 
the intersection of SU(4) x SU(2) x U(l) and SU(5) x U(l) in SU(6), which is the 
Standard Model gauge group apart from an extra U(l). Matter fields arise from 
5d bulk fields. This model is sketched in Fig. [2| 

Regarding SU(6) as a spontaneously broken global symmetry of the 4d the- 
ory, the Higgs fields are massless in the absence of SUSY breaking because 
their imaginary parts are Goldstone bosons associated with the breaking to 
SU(4) X SU(2) X U(l), and their real parts are then protected by SUSY. SU(6) is 
also broken explicitly by boundary conditions on the Planck brane (or on the level 



In 13 this chiral adjoint is denoted by S. 
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UV brane 



$ 



IR brane 



Figure 2: A warped holographic GUT model : $ (light green dashed curve) is an IR 
brane field. Matter fields are localized towards the UV brane, only slightly so for the 
third generation but more pronounced for the first two generations. 



of the 4d theory by gauging only its SU(5) x U(l) subgroup), and therefore the 
Higgs fields are merely pseudo-Goldstone bosons. However, this explicit break- 
ing gives tree-level masses to only 12 out of the 16 Goldstone modes. In other 
words, if SU(5) x U(l) C SU(6) is gauged, only 12 of the Goldstone bosons will 
be eaten by the Higgs mechanism, corresponding to the 12 broken generators in 
SU(5) X U(l) -^ SU(3) X SU(2) x U(l) x U(l). The remaining four Goldstone 
modes and their complex partners form two massless weak doublets which are 
identified with the MSSM Higgs fields. (For some earlier work along these lines, 
see |T,^.) 

Provided that the SUSY breaking mechanism respects this symmetry breaking 
pattern, it will lead to Eq. ^ because once more the fields from the combination 
$''' + $ can pick up a tree-level mass term, while the (pseudo-) Goldstone bosons in 
$^ — $ remain massless. Since the SU(6) symmetry is explicitly broken, radiative 
corrections can again lift the relations of Eq. ^ below the SUSY breaking scale, 
as required for phenomenology. 

Localizing the Higgs fields on one of the branes rather than in the bulk natu- 
rally allows for the possibility that they are "sequestered" from the supersymmetry 
breaking fields Z\ which could be confined to the other brane. Writing the ef- 



fective quadratic Higgs Lagrangian as in Eq. (16), this means that there exists a 
frame in which 

gF' = 0. (18) 

In that case any tree-level Higgs mass terms should arise purely from gravitational 
effects. We are parametrising these by the chiral compensator (f. Note that this 
type of scenario includes the case of radion mediation, i.e. taking the dominant 
source for SUSY breaking to be the F-term of the radion modulus, since one can 
always redefine ip such that the radion does not couple to IR brane fields. It is well- 
known that F^ alone cannot induce scalar soft masses classically (in accord with 



Eq. (17)), with the leading contribution to mj^^ and mj^^ coming from anomaly 
mediation and thus suppressed by a loop factor p!4, 15 . F'^ will however give rise 
to a /x and a S^ term at tree-level, 
and (flTl): 

J- 1 



l/^l 



l^^l 



\B,, 



One thus obtains a special case of Eqns. ^ 

(19) 



\m 



H^ 



< IF^I 



Moreover, in the holographic GUT model with SUSY breaking on the UV brane 
there is no reason to expect the matter soft terms for the first two generations 
to be suppressed. The wave functions of the first and second generation have in 
fact similar to slightly larger overlaps with the UV brane, as compared with those 
of the third generation. This is because we are generating the Yukawa hierarchy 
from hierarchical wave function values on the IR brane, as depicted in Fig. [2} 

Needless to say, this model also allows for supersymmetry breaking on the 
IR brane. In that case Higgs soft masses will generically be induced by direct 
contact interactions, and the ordering of sfermion masses will roughly correspond 
to the Yukawa hierarchy. In fact in general one can have a variety of different 
contributions to SUSY breaking mediation from brane-localized fields, as well as 
from the radion and from Lp. 



3 Soft term patterns 

In the common public SUSY spectrum generators, /i and 5^ at the GUT scale are 
computed according to their renormalization group equations from their IR-scale 
values, which in turn are calculated from rriz and the given tan^S. Here we use 
SOFTSUSY [16] . Following [4|, we implement the Higgs mass relation Eq. ^ by 
iteratively adjusting the Higgs soft masses, 

^^1 = ^k -^ ^hB^ - |/i|^ (20) 

at Mx = Mgut- Here e^f = ±1 takes care of the sign ambiguity in 5^ which we 
mentioned in the previous section. Throughout this analysis we assume gaugino 
mass unification. Mi = M2 = M3 = M1/2 at Mqut- The free parameters in our 
study are thus Mi/2(Mgut), tan/3(M^), the two signs en and sign(/x), and the 
sfermion soft terms at Mqut- We take the latter to be fiavour- diagonal. 

Whether EWSB and a viable phenomenology can be obtained strongly depends 
on the sfermion soft terms. Two limiting cases are of particular interest: 

• assuming a common sfermion mass mo and a common trilinear coupling 
Aq. This makes the DHMM models a subclass of non-universal Higgs mass 
models (see e.g. [17 and references therein) with m|^^ = m^^^ ("NUHMl" in 
the terminology of Ts]); 

• no-scale boundary conditions for the first and second generation, mo(l, 2) = 
Ao(l, 2) = 0, but allowing for arbitrary soft terms in the third generation. 

The first case may be considered as representative for a generic scenario with 
all sfermion soft terms of the same order of magnitude, whereas the second case 
represents models with hierarchical soft terms refiecting the Yukawa hierarchy. 
As we have seen in Section [2} both these scenarios are well motivated from the 
model-building point of view. Moreover, it is interesting to investigate whether 



the stronger condition of Eq. (19), m|^^ ^ -^ and |/x| -^ |S^|, can be realized. 



It turns out that the following patterns emerge in the soft terms: 

8 



1. in almost all of the admissible regions of parameter space, ch corresponds to 
the GUT-scale sign of 5^; 

2. for sign(/x) = +1, B^ has almost always the same sign as At at the GUT 
scale (and the opposite one if sign(/x) = — 1); 

3. for sign(/i) = +1, the stricter relation 

^ffi = ^k = 0' ^H B^ = |/i|^ at Mgut (21) 

can only be satisfied with 6/f = +1. 

These observations can be explained by a close inspection of the relevant RGEs, 
as we will now detail. 



To start with, it is useful to recall the dominant contributions to the one-loop 
RG evolution of the stop trilinear At. These involve At itself and the gluino mass: 

>±^. .., ,. . . 32 

The large gluino contribution will drive At to a large negative value towards "late 
times" (i.e. towards the low energy scale), until it is compensated for by the first 



d 32 

167r2-^ = 12 Ivtl^At + -gJMs + ... (22) 



term in Eq. (22) 



Now concerning point 1., if sign(5^) does not match with ch at the GUT 
scale, this implies that the GUT-scale ml and ml are negative. The running of m^ 
and 7712 i^ almost exclusively due to the running of the soft masses m|^^ and m^^^ ^ 
since /x is approximately constant. The dominant terms in the one-loop RGE for 
the up-type Higgs soft mass-squared are 

IGTr^^mf,^ = 6 \yt\^ {\A\^ + m\^ + m%^ + ml^) - ^gl\M^f -r . . . . (23) 

where m^^ and m?Q^ are the soft masses of the third generation up- type squarks and 
squark doublets respectively. In scenarios like the CMSSM, by scalar mass univer- 
sality the terms in parentheses are typically positive. Thus mt^ is driven to lower 
values as the RG scale decreases, assisted also by the top Yukawa coupling and 1^4^ | 
growing large. Eventually radiative electroweak symmetry breaking is triggered. 
Most of the DHMM parameter space also has this property. There is only a tiny 
region with initially small A^, very large negative m|^^ ^^ small or negative squark 
soft masses-squared, and sizeable M1/2, in which mt^ runs up significantly at first. 
In that case it can be driven to positive values even if it starts out negative, and 
electroweak symmetry breaking can be triggered later when the At contribution 



in Eq. (23) dominates and when also the squark masses have grown positive. 



Concerning point 2. above, note that By, at the low scale should be somewhat 
small compared to the typical soft masses. This is in order to satisfy the elec- 
troweak symmetry breaking conditions Eq. ([2]), and in particular to have at least 
moderately large tan/5 (tan/3 > 5 say). Let us for now assume positive /x. The 
one-loop RGE of By, is dominated by the At and gaugino contributions: 

167r2-5^ = 6/x A \yt\^ + QiiM2gl + ... (24) 



The gaugino contribution tends to dominate the RG evolution of 5^ at scales close 
to the GUT scale, driving 5^ down. However, eventually At itself will run large 



and negative because of the gluino contribution to Eq. (22). Far in the IR it will 
thus primarily drive the S^ evolution, causing Bj^ to run up instead. For sizeable 
and positive initial GUT-scale At^ this latter effect will be less important. But 
for small or even negative GUT-scale values. At will quickly evolve towards large 



negative values, and thus dominate over the gaugino term in Eq. (24). 



We conclude that positive starting values for the GUT-scale At are preferred if 
B^> Q (in which case B^ should mainly run down towards the electroweak scale, 
to end up small) and small or negative values are preferred if S^ < (in which 
case it should mainly run up, to end up positive). For negative /x, the signs are 
reversed. These correlations are illustrated in Fig. [3] with the help of some slices 
through the parameter space, at fixed sfermion and gaugino masses, fixed tan/3, 
and universal trilinears. 

Our choice for the soft masses in Fig. |3] may seem rather high, but as can be 
seen from Fig. [4l satisfying the LEP Higgs bound requires fairly large M1/2, the 
more so the larger A^ is. 

Finally let us come to point 3.: The At parameter also enters the ttz^^^ RGE 
Eq. (23). A large 1^4^ | will accelerate the decrease of Tri\j^ when running down 



from the GUT scale, which is of course a particularly severe effect if At starts out 
negative, i.e. if en = — sign(/i). Then m^^^ ^iU ^^^ negative too quickly, unless 
there is some other contribution to counterbalance the effect of At. Large gaugino 
masses could provide such a contribution, but they would again slow down the 



S^ evolution, as is evident from Eq. (24). If we allow for non- vanishing Higgs 



soft masses-squared, they can in particular be negative and thus counteract the 
At effect in Eq. (23). However, in the case that ^1/^2 ^^ constrained to vanish. At 



should be positive at the GUT scale to minimize its effect on the Tri\j^ running. In 
addition, small (or even negative) GUT-scale squark masses-squared are preferred 
to slow down the m^^ evolution. Indeed, in Figs, p^ and the case m|^^ ^ ~^ 
occurs only for en = sign(/i) and requires large positive A^ to start with, leading 
to small negative At at the EW scale. 

As can also be seen from Fig. |4| large positive ^40 leads to a tension with 
the direct search bound from LEP of rrih > 114: A GeV at 95% C.L. |19|, even 
when taking into account a 2-3 GeV theoretical uncertainty [20] on the calculation 
of rrih in the MSSM. This can be understood as follows: At lowest order, the 
light CP-even Higgs boson of the MSSM is at most as heavy as the Z^ boson, 
^^ ^ ^1 cos^ 2/3. Radiative corrections have to lift rrih above the LEP limit. The 
dominant effect is proportional to the fourth power of the top Yukawa coupling, 7/^, 
and comes from an incomplete cancellation of top and stop loops. This increases 
rrih approximately to 



Sg^m^ 



"^^^"^^-8^^ 



w 



2\ V2 / V2 ^ 

(25) 



miMU^r,-^^ 



m^ I ^'^^ 



Ml \ 12M|y_ 



where 



M| = J(m|+m|), Xt = At-ficot(]. (26) 
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Figure 3: A set of slices through the four branches of DHMM parameter space with 



M 



1/2 



1 TeV, tan/3 = 10, and ttiq = 500 GeV. For the panels on the left, regions 



of larger Aq do not lead to electroweak symmetry breaking, as can be seen from the 
pseudoscalar Higgs mass tti^ approaching zero. The same is true for regions of smaller 
Ao for the panels on the right. The upper panels show, as explained in the text, that 
for sign(/x) = +1 the GUT-scale sign of At is equal to the GUT-scale sign of B^ in 
almost all of the allowed regions. The lower panels show that this correlation is reversed 
if sign(/x) = —1. As also explained in the text, 7n|^^ ^ = is only possible if e/f = +1 
and sign(/x) = +1 (top right panel). Finally note that there is a tiny slice of the allowed 
parameter space where e^f does not correspond to sign(S^) at the GUT scale. In this 
region l^tl is small and ttt.^^^ is large and negative as expected. 
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Figure 4: Contours of constant m^^ ^ (full red lines) and |/x| (dashed blue lines) at Mqut 
in the M1/2 versus Aq plane, for tan^ = 10 and mo = 500 GeV. Grey regions do not lead 
to electroweak symmetry breaking, light yellow regions are excluded by 6 ^ 57 (at 2cr), 
and dark yellow regions have rrih < 114 GeV, with the lines indicating rrih = 111 GeV. 

For large tan^S and large |/x|, also bottom and sbottom loops become important, 
giving an analogous contribution proportional to y^. For details see, e.g., [21] 
and references therein. The logarithmic sensitivity to the average stop mass Ms 



in Eq. (25) suggests that heavy stops are preferred in order to render rrih large 



enough. However, this sensitivity is rather mild, and the dependence on the stop 
mixing parameter Xt can be at least as important. Indeed, rrih initially increases 
with \Xt\ and reaches maximal values for Xt = ±VQMs] this is known as the 
'maximal mixing' or m^^^ case, see again |21|. Therefore a large low-scale \At\^ 
together with moderately large tan/3, is favoured to satisfy the LEP Higgs mass 
bound. This is exactly what we find in the right-hand side panels of Fig. [4| For 
too large starting values of At^ the low-scale \At\ will be too small (recall that 
At generically runs towards negative values) and the Higgs mass bound becomes 
important. 
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Figure 5: Higgs mass parameters in the MSSM with non-universal Higgs masses, for 
a universal gaugino mass M1/2 = 1 TeV, a universal sfermion mass mo = 500 GeV, 
tan/3 = 10, and vanishing Higgs soft masses. Note that there is only a "DHMM point" 
(with meeting |/xp and S^ curves) for /i > 0. 



The DHMM model with universal sfermion soft terms is a special case of a 
NUHM model. Fig. [5] shows for comparison the dependence of m^, /x^ and 5^ 
in a general NUHMl model with mj^^^ = 0. The other parameters are as in 
Fig. 3I The CMSSM limit with m^^ ^ = ml gives almost the same picture, the 
only difference being a slightly larger m\ and slightly smaller /i^. Note that there 
is only one "DHMM point" in Fig. p^ (with meeting |/xp and 5^ curves), which 
occurs for /x > 0. Away from this point where the models coincide, DHMM has 
a much larger |/i| and smaller ttia than NUHMl (or the CMSSM), cf. Fig. |3| In 
particular, in the DHMM case tua becomes small for small 1^40 1, and we can have 
rriA ^ Ml/2 even for small tan^. This will be important later when we consider 
the neutralino relic density. At this stage we just remark that in Fig. [4]5-channel 
annihilation through the Higgs funnel occurs for small |ylo| ^ 100 — 300 GeV. 

Let us now explore the consequences of non-universal third generation soft 
terms. This is interesting in particular for the gauge-Higgs unification models 



discussed in Section |2.H where we expect vanishing first and second generation 
soft masses, mo(l,2) ^ 0. In this case we have roughly m|^ ^ (0.39 Mi/2)^ — 



0.052 S'gut, and mj^ ^ (0.68 Mi/2)^ + 0.026 S'gut, which has to be compared 
to m^o ^ 0.43 Ml/2. Here, S'gut is the GUT-scale value of the hypercharge S 
parameter. 



S 



{ml^-ml^] 



+ Tr(mg — 2mfj + m|) 



+ ml- ml) 



(27) 



We see that a non-zero and negative Squt of about — (0.8Mi/2)^ to — (3.3Mi/2)^ 
is necessary if one wants the neutralino to be the lightest SUSY particle (LSP). 

in DHMM, the way to ensure a neutralino LSP is 



Since we have m|^^ 



m 



H2 



non-universality of the third generation, as illustrated in Fig. 6][^ We can see 



^We use the notation mu^ 



m 



Us 



'4\ 



mj 



so the sign of mus is actually that of m'^ 



Us 



and 
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Figure 6: Contours of constant m^^ ^ (full red lines) and \/j.\ (dashed blue lines) at Mqut 
in the TTij/g versus ruQr^ plane, for M1/2 — 1 TeV, tan/3 = 10 and vanishing lst/2nd 
generation soft terms. Moreover, /x > 0, e/f = ±1, and Aq = zbMi/2 (upper row) and 
Aq = zh3Mi/2 (lower row). Back regions do not lead to electroweak symmetry breaking, 
gray regions have a slepton LSP, and white regions a neutralino LSP; the light yellow 
stripe is excluded by 6 ^ 57 at 2cr, while the dark yellow stripes have ruh < 114 GeV. 
/i < gives qualitatively very similar results. 



that taking a slice along mu^ = ttiq^ = tuds^ as we have done in the previous 
plots, indeed results in qualitatively similar patterns as the more general case in 
which the squark soft masses are split. On the other hand, negative soft masses- 
squared can give much smaller mj^^^ and |/i|. This is of interest in particular for 
^H = sign(/x) = +1, where one can achieve a mixed bino-higgsino LSP (see the 
dark matter discussion in the next section) . Note moreover that in the RHS plots 
of Fig. pi mjj^ 2 is very sensitive to Aq^ while /x does not vary much when passing 



analogously for tuq^ etc. 
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from Aq = Ml/2 to ^40 = 3Mi/2. This is in accord with Fig. pj 

We conclude this section with a few remarks on potentiaUy dangerous tachyonic 
directions. In our analysis we have permitted tachyonic GUT-scale masses for both 
the Higgs and the sfermion fields. This is well-known to generally lead to charge- 
and colour-breaking minima in the potential, as well as to directions in field space 
which are unbounded from below (at tree-level and without higher-dimensional 
operators); see, for instance, |22|. 

While tachyonic scalar masses appear to rule out a large part of the parameter 
space at first sight, two points must be considered. First, a careful analysis is 
required in each case to determine if such dangerous vacua really are present. 
In particular, calculations using only the RG-improved tree-level potential may 
give unreliable results if the field VEVs are found to be vastly different from 
the renormalization scale, because of the presence of large logarithms in the loop 
corrections. Second, these vacua need not be dangerous even if they are present. 
If the tunnelling rate from our false vacuum is sufficiently small, our vacuum may 
well be effectively stable on cosmological timescales. It then depends on early- 
universe cosmology whether or not it is preferred for our universe. See i23i for 
a recent analysis of the CMSSM and of NUHM models in that context, and (24] 
for a recent analysis of the cosmological lifetime of related Higgs-exempt no-scale 
models. 

A detailed investigation of charge- and colour-breaking minima is beyond the 
scope of this work. We therefore merely stress that we expect them to appear 
in large regions of parameter space, but depending on their lifetime and on the 
cosmological scenario, these regions may still be acceptable phenomenologically. 



4 Markov Chain Monte Carlo analysis 

So far we have only considered the constraints from rrih and b -^ 57, and taken 
one- or two-dimensional slices through the parameter space. In order to take into 
account more constraints and in particular to find regions of parameter space where 
the neutralino LSP is a good dark matter candidate, we next perform a Markov 
Chain Monte Carlo scan of DHMM models. As above, we consider the two cases of 
(i) universal sfermion soft terms and (ii) vanishing first/second but non-universal 
third generation soft terms. 

MCMC is an efficient method to probe a large-dimensional parameter space, 
and to gain information about it by using Bayesian statistics. The basic idea 
is to set up a random walk, starting at some parameter point and proposing a 
candidate next point at random nearby. This candidate point is then accepted 
or rejected at random, with an acceptance probability depending on its likelihood 
compared to the likelihood of the original point. Parameter points which are more 
likely to reproduce existing experimental data and constraints within errors have 
a greater probability of being accepted. If accepted, the new point is chosen as 
the starting point and the procedure is iterated. Otherwise it is repeated with the 
old starting point. A properly set-up ensemble of Markov chains should eventually 
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fill out all the allowed parameter space, with a high density of points in those 
regions which are best compatible with existing measurements. In the sense of 
Bayesian statistics, the distributions of points are interpreted as probability density 
functions. MCMC provides a simple means to marginalise these distributions and 
to evaluate probability regions. 

The setup and procedure of our MCMC analysis closely follows [25 1, and we 



refer the reader to this paper for technical details (see also 26 -[29]). Here we 
just explain the constraints and priors used in our analysis. We apply the limits 
from direct SUSY |30| and Higgs 19,31] searches at LEP. The computation of rrih 



suflFers from a theoretical uncertainty which has been estimated to amount to up 
to 2-3 GeV [20J. This theoretical error is most likely non-Gaussian and can give 
an underestimation as well an overestimation of rrih- We therefore use the direct 
experimental search limit for a SM-like Higgs of rrih > 114.4 GeV at 95% C.L. 
without further modification. One should however bear in mind that the favoured 
regions of parameter space may in fact be somewhat larger (where the Higgs mass 
is underestimated by the calculation) or smaller (where it is overestimated) than 
the ones we find. 

Regarding the anomalous magnetic moment of the muon, we limit our scans 
to // > 0, which gives a positive SUSY contribution, but do not require that SUSY 
explains the discrepancy between the measurement and SM prediction; instead we 
only apply an upper limit on Aa^^^ . 

The complete set of constraints applied is given in Table [T] For observables on 
which there is merely an experimental upper or lower bound available, we use a 
Fermi likelihood function Li. For quantities which have been measured, we use a 
Gaussian likelihood function L2. The total likelihood of a parameter point is the 
product of all individual likelihoods, L = nn(^0^- -'-^ ^^^ notation of Table fll we 
have 



Li{x,Xo,dx) = -— vm ^ L2{x,Xo,dx) = exp 

1 + exp[(x — Xo)/dx\ 



(x - Xpf 
2dx^ 



. (28) 



The neutralino relic density, the B-decay branching ratios, Aa^^^^, and the SUSY 
mass limits are evaluated with micrOMEGAs 39,40|j^ 

We choose to work with two different prior probability distributions. Our first 
prior is fiat in the GUT-scale soft parameters and in tan/3. That is, within a 
certain fixed range, any value for a given parameter is treated as equally probable. 
As a second prior, for comparison, we use a "naturalness prior" |26j: Since a prior 
choice ultimately refiects theoretical prejudice as to what parameter choices should 
be more or less likely, we find it appropriate to use a prior which disfavours the 
more fine-tuned parameter points. The main source for fine-tuning in the MSSM 
is caused by the sensitivity of the elect roweak scale to parameter variations. We 



^In the likelihood function for Qh^, we use the 2008 central value of 38 with a Gaussian width 
of about 10%. This is to approximately account for uncertainties from the cosmological model, 
from the data sets used, and from the SUSY spectrum calculation. It is consistent with the most 
recent determination of Qh^ from seven-year WMAP data, published in early 2010 l4ll|42|. 
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Observable 


Limit 


Likelihood function 


Ref. 


ruh 


> 114.4 


Li(a;, 114.5, -0.6) 




31 




rrit 


173.1 ±1.3 


L2(x, 173.1, 1.3) 




32 




mw 


80.398 ± 0.025 


L2(x, 80.398, 0.025) 




te 




BR(6 -^ 57) 


(3.52 ± 0.34) X 10-4 


L2(x,3.52x 10-^0.34x 10-^) 




34 


,35 




BR{Bs -^ /i+yU") 


< 5.8 X 10-8 


Li(a;,5.8x 10-^5.8x 10-^°) 




36 




R(5„ ^ TUr) 


1.11 ±0.52 


L2(a;, 1.11, 0.52) 




34 




AafSY 


< 4.48 X lO-'^ 


Li(a;,4.48x 10-^4.5x 10-^^) 




37 




nh^ 


0.1131 ±0.0034 


L2(x, 0.113, 0.011) 




M 




SUSY mass limits 


LEP limits 


1 or 10-*^ 




ffl 





Table 1: Observables used in the hkelihood calculation. Li and L2 are defined in 



Eq. (28) 



therefore use a fine-tuning measure c defined as |43| 



max 



dlnmz 



d In ttj 



(29) 



Here {a^} includes all GUT-scale soft masses and trilinear soft terms, as well as /x. 
With the naturalness prior, every parameter point is then weighted with a measure 
1/c, thus penalizing the more fine-tuned ones. 

Before presenting the results, let us comment on the viable dark matter regions. 
In general in the MSSM there are only a few mechanisms that provide the correct 
amount of neutralino annihilation consistent with cosmological observations (see 
e.g. |44| for a review). In the DHMM case we consider here, we expect: 



A) Coannihilation with sleptons. This requires small neutralino-slepton mass 
differences of roughly 10 — 1 GeV for m^o ^ 100 — 500 GeV; for heavier 
LSPs, coannihilation with sleptons alone is not efficient enough. Another 
possibility is coannihilation with light ti or 6i, which is efficient for larger 
mass differences, or larger LSP masses. 



B) Annihilation through 5-channel pseudoscalar Higgs exchange. Here the key 
quantities are the distance from the A pole, ttia — 2m^o, and the width of the 
A resonance. The process is efficient for a bino LSP, although some higgsino 
admixture is necessary to provide the XiXM coupling. 



C) Annihilation of a mixed bino-higgsino LSP through t-channel chargino and 
neutralino exchange, and through 5-channel Z exchange. This requires a 
sizable LSP higgsino fraction fn ^ 0.25%. Heavier LSPs need a larger 
higgsino fraction, so that eventually coannihilation with other neutralinos 
and charginos also becomes important. Besides, if 2m^o ^ m^, 5-channel A 
exchange also contributes in this region. 



^xi 



Finally note that throughout this work the squark and slepton mass matrices 
and 74-terms are assumed to be diagonal. The issue of ffavour-changing neutral 
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Figure 7: Marginalized ID posterior probability distributions of the input parameters 
for universal soft terms and e/f = +1. Black lines are for flat prior, red lines for natural 
prior. Dimensionful quantities are in GeV. 
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Figure 8: Same as Fig. [t] but for en = —1; black (blue) lines are for flat (natural) prior. 



currents due to non-diagonal terms arising in particular in the warped case [45j-[47 
is left for a separate work |48|. 



4.1 Results for universal soft terms 

Here the model parameters to scan over are universal gaugino and sfermion mass 
parameters M1/2 and mo, a universal trilinear coupling Aq^ and tan^S. In addition 
there are the two discrete parameters sign(/i) and ch- We choose /x > as favoured 
by BR(6 -^ 57) and run ten chains with 10^ iterations each, for both e/f = +1 
and —1, allowing M1/2 to vary from to 2 TeV, ttiq from to 5 TeV, Aq within 
±10 TeV and tan/3 from 2 to 60. 

Figure [7| shows the marginalized ID posterior probability distributions of the 
input parameters comparing flat (in black) to natural (in red) prior. The case of 
6/f = — 1 is shown in Fig. [S] As can be seen, in both cases the naturalness prior 
results in a pull towards smaller masses and smaller tan^. The general features, 
which are detailed below, however remain the same. 

Ml/2 is bounded from below by the Higgs and SUSY mass limits, and from 
above by the requirement of suflicient neutralino annihilation. The processes that 
bring the neutralino relic density within the desired range are A) or B) from above: 
coannihilation with sleptons {cr^ Jir^ or fi) or annihilation through the Higgs 
funnel. On the other hand, we do not find any region where the LSP higgsino 
fraction is large enough to render processes C) eflicient. Coannihilation with stops 
or sbottoms is also absent. For e^f = +1 it becomes diflicult to achieve small 
enough |m^ — 2m^o| and m^~ — m^o when m^o > 750 — 800 GeV. For ch = —1 this 
is the case when m^o > 600 GeV. 
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Figure 9: Contours of 68% and 95% probability in the BR(6 -^ 57) versus tan/3 plane, 
on the left for en = +1, on the right for en = —1- The green shading maps the average 
likelihood per bin, normalized to the maximum likelihood. 

The relic density constraint also prefers higher tan/3, for which the Higgs 
funnel is more efficient. This is the reason for the preference of high tan/3 in the 
distribution for 6/f = +1 and flat prior (which is still softened by the natural prior). 
High values around tan/3 ^ 50 are constrained by BR(6 -^ 57) becoming too low. 
For e/f = — 1, the tan^ distribution is more flat because the BR(6 -^ 57) constraint 
becomes effective earlier as tan f3 grows. These correlations are illustrated in Fig. |9) 

Regarding the mo probability distribution, the peak at low mo is where coan- 
nihilation with sleptons takes place. Slepton coannihilation, and with it the low 
mo peak, becomes more relevant when using the natural prior because of its prefer- 
ence for smaller tan^S for which the Higgs funnel is less efficient. Overall, however, 
the annihilation through the pseudoscalar resonance is by far the dominant mech- 
anism: for 6/f = +1 and flat (natural) prior, 88% (83%) of the points exhibit 
predominantly annihilation into 66, while 10% (15%) predominantly show coan- 
nihilation with sleptons. For 6/f = — 1, the A resonance is more diflicult to hit, 
partly because tan^S is smaller, so that for flat (natural) prior 22% (26%) of the 
points predominantly show slepton coannihilation. 

As opposed to the CMSSM there is no "focus point" behaviour in this scenario: 
The mo distribution shows a clear preference for lower values < 2 TeV, and mo is 
signiflcantly correlated with /i. In fact the CMSSM focus point hinges on having a 
single parameter which governs the scalar soft masses for both Higgs and matter 
fields. This is clearly not the case in DHMM models. 

Finally, the Aq distribution confirms our discussion of the sign correlations in 
Section [3l 

An important issue in our considerations are the values of mj^^ ^ , fi and 5^ 
at the GUT scale resulting from the DHMM condition. 



we 



In Figs. [To] and 11 

therefore show 2D posterior probability distributions for these parameters jjThe 
tight correlation between mj^^ ^ and // is clearly visible. Moreover, as can be seen. 



for both 6/f = ±1 small 



m 



^1,2 



and |/x| prefers small values of M1/2, tuq and chAq. 



Regarding consequences for experiments. Figs. [12] and [13] show ID posterior 
probability distributions for SUSY and Higgs masses. Also shown are the distribu- 



^To limit the proliferation of figures we only show 2D distributions for flat prior; those for 
natural prior look very similar. 
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Figure 10: Contours of 68% and 95% probability showing correlations between ^n^^^^, 
/i, Bjj, and the input parameters for universal soft terms, en — +1 and flat prior. 
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Figure 11: Same as Fig. 10 but for e/f = — 1. 



tions for the LSP higgsino fraction fn and the cross section for spin-independent 
direct detection a^^. The pull of the natural prior towards lighter masses and 
in particular towards smaller /x is again evident. We also note that most of the 
parameter space lies within reach of the LHC at 14 TeV centre-of-mass energy. In 
fact, for ch = +1 (—1) and flat prior, 82% (97%) of the points have gluino and 
squark masses below 3 TeV. Moreover, 55% (58%) of these points have sleptons 
that are lighter than the X25 so that a same- flavour opposite-sign dilepton signal 
from X2 ~^ i^i^ -^ ^^^^Xi ^^y be visible in SUSY cascade decays (if decays 
into sleptons are absent or kinematically suppressed, then X2 ~^ ^Xi is the most 
important decay mode of X2)- F^^ naturalness prior, 88% (99%) of the e^f = — 1 
points have gluino and squark masses below 3 TeV, with 64% of these featuring 



rrie^f < rn^o. 



If the X2 decay into sleptons is open, X2 ~^ e^e^,/x^/x^ has up to about 
40% branching ratio. It is however important to keep in mind that owing to the 
universality assumption, the typical mass ordering is rrif^ < rriej^ < rrie^- Therefore 
X2 -^ f^T^ decays are often dominant. 
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Figure 12: Posterior probability distributions of the most relevant masses for universal 
soft terms and ch — +1- The bottom-right plots show the LSP higgsino fraction, 
Ih •= |^i3p + I^i4p5 and the spin-independent scattering cross section on protons. As 
above, black lines are for flat and red lines for natural prior. 
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Figure 13: Same as Fig. 12 but for en = —1; black (blue) lines are for flat (natural) 
prior. 
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Figure 14: 2D posterior probability distributions of the fine-tuning measure c for natural 
prior. The contours enclose regions of 68% and 95% probability. The top row (red 
contours) is for e/f = +1, the bottom row (blue contours) for e/f = — 1. 



Concerning direct dark matter detection, we note that because the LSP is 
always almost a pure bino, the neutralino scattering cross section on proton is 
typically of the order of 10~^^ — 10~^° pb and hence beyond the reach of current 
experiments. 

Finally, we observe that even with the natural prior the fine-tuning tends to be 
very large, of the level of per-mil, and points with c < 100, corresponding to less 
than 1% fine-tuning are difficult to obtain. Correlations of the finetuning measure 



c are illustrated in Fig. [M] for natural prior. The lowest fine-tuning occurs for 
small Ml/2, medium tan/5 ^ 20-30, ^40 ^ and mo ^ 1 TeV, with // being around 
1.5-2 TeV. 



4.2 Results for vanishing lst/2nd generation soft terms 

Let us now turn to the pattern of soft terms obtained from models such as the 



gauge-Higgs unification (GHU) model of Section 2.1 Here the first- and second- 
generation matter fields were localized on a brane and SUSY breaking was me- 
diated by the radion, leading to vanishing lst/2nd generation and non-universal 
3rd generation soft terms. We will call this "GHU-like boundary conditions" in 
the following. The free parameters in this case are M1/2, tan^S, and the third- 
generation soft terms mgg , m^/g, vfiB^^ At^ Aij^ ^Lg, '^£;3, A^. We allow M1/2 to 
vary from to 2 TeV, tan/3 from 2 to 60, '^gg, [73,^3 within ±25 TeV^, Tn\^^^^ from 
to 4 TeV^ and A,6,r within ±10 TeV. 

The marginalized ID posterior probability distributions of the input parame- 
ters are displayed in Fig.fT5]for en = +1 and in Fig. [TBlfor en = —1. Analogously, 



Figs. [17] and [18] show the probability distributions of masses, jjl parameter, LSP 
higgsino fraction, and the spin-independent LSP scattering cross section on pro- 
tons. 
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Figure 15: Marginalized ID posterior probability distributions of the input parameters 
for GHU-like boundary conditions and e/f = +1. Black lines are for flat prior, red lines 
for natural prior. Dimensionful quantities are in GeV. 
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Figure 16: Same as Fig. 15 but for en = —1; black (blue) lines are for flat (natural) 
prior. 



Two important differences to the case of universal soft terms are that M1/2 
can now go to much higher values, and that tan^S peaks around 10. The reason 
is on the one hand that due to the no-scale boundary conditions for the lst/2nd 
generation, coannihilation with selectrons and smuons becomes more likely; this 

< 500 GeV. Accordingly, there are distinct peaks at 
and 



is mainly relevant for m^o 



m^o ^ 400 GeV in Figs. 



17 
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corresponding to the peaks at M1/2 ^ 900 GeV 
in Figs. [15] and [T6| On tKe^othenTand, due to the non-universal 3rd generation we 
can obtain smaller values of /x, and hence processes C) become important. This 
is mainly relevant for heavy Xi ^^d leads to the peak at large M1/2 for 6/f = +1. 
For ch = —i^ 1^ tends to be larger {i.e. fn tends to be smaller) and consequently 
the high M1/2 region is less favoured. Besides, for both e^f = ±1, we find some 
coannihilation with bi and/or ti, though this is diminished by the naturalness 
prior. (For ch = — 1, this leads to the peak at large negative ttidsi which gives 



light 61 ^ 6i^, c.f. Fig. [18| Coannihilation with ti is less frequent, in particular for 
e/f = +1, as the Higgs mass bound pushes the stop masses up.) 

Concerning collider phenomenology, we first observe that, because of the van- 
ishing lst/2nd generation soft terms, the X2 ~^ £^£^ -^ ^^^^Xi decay, with £ = e 
or /i, is almost always present. Staus are heavier and hence much less important 
for the X2 decays. Second, for e^f = — 1, squark and gluino masses peak around 
2 TeV, which means that the LHC at 14 TeV centre-of-mass energy has again a 
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very good discovery potential over most of the parameter space. More precisely, 
88% of the 6if = — 1 points have rriq^g < 3 TeV. For ch = +1^ on the other hand, 
we find that a considerable fraction of the parameter space lies beyond the reach of 
the LHC. In this region the Xi is heavy and is very likely to have a large higgsino 
fraction (since we require Qh'^ ^ 0.1). In turn this leads to a large cross section 
for direct dark matter detection of up to around 10~^ pb, see the bottom right 



plot in Fig. [T7| Interestingly, this is just at the edge of current CDMS-II exclusion 
limit [49 1 for heavy masses H The 2D probability distributions in the plane a^^ 
versus m^o are shown in Fig. fT9| for the natural prior. It is very gratifying that 
these models can be experimentally tested with complimentary methods, by both 
LHC and direct dark matter searches. 



For completeness. Figs. 20 to [23] show various parameter correlations in 2D. 
It is interesting to see that mj^^ ^{Mqut) = is easily obtained for ch = +1, but 
does not occur for 6/f = — 1. Moreover, the sign correlation between e^f, 5^ and 
At discussed in Section |3] is evident. 



^ While we have not used constraints from direct dark matter searches in the MCMC, a 
posteriori it turns out that only about 1% of the points with higgsino LSP violate the current 
limits. 

24 




400 600 

m 





1000 2000 3000 4000 5000 

m + 







1000 2000 3000 4000 5000 



1000 2000 3000 4000 5000 



1000 2000 3000 4000 5000 



1000 2000 3000 4000 5000 






|°giofH 




Figure 17: Posterior probability distributions of the most relevant masses, //(Mew)? LSP 
higgsino fraction, and the spin-independent LSP scattering cross section on protons for 
GHU-like boundary conditions with e/f = +1. As above, black lines are for flat and red 
lines for natural prior. 
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Figure 18: Same as Fig. 17 but for en — —1; black (blue) lines are for flat (natural) 
prior. 
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Figure 19: Probability distributions in the plane a^^ versus m~Q for naturalness prior, 
on the left for e^f = +1, on the right for e^f = — 1. The inner (outer) contours enclose 
regions of 68% (95%) probability, the green shading maps the average likelihood, and 
the black lines show the limit from CDMS-II, which is currently providing the strongest 
bound in this mass range. 
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Figure 20: Contours of 68% and 95% probability showing correlations between the most 
relevant input parameters for GHU-like boundary conditions, ch — +1 and naturalness 
prior. 
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Figure 21: Same as Fig. 20 but for e/f = — 1. 
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Figure 22: Contours of 68% and 95% probability showing correlations between rnj^^ ^, /i, 
Bj^ and the most relevant input parameters for GHU-like boundary conditions, e/f = +1 
and naturalness prior. 
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Figure 23: Same as Fig. 22 but for e/f = — 1 
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5 Conclusions 

Among the many possible embeddings of the MSSM into a grand-unified theory, 
there are some interesting classes of models which predict a degenerate GUT-scale 
Higgs mass matrix. We have investigated the origin of this prediction in some 
example high-scale models, as well as its consequences for low-scale mass spectra 
and phenomenology. 

With the additional assumption of universal GUT-scale gaugino masses (which 
is valid in most simple GUT scenarios) the low-energy spectrum still depends 
sensitively on the sfermion soft terms. Different high-scale models will give rise 
to various patterns of sfermion masses and trilinear terms. We have chosen to 
investigate two representative cases in detail: first, universal sfermion soft terms, 
and second, vanishing soft terms for the first two generations but non- vanishing 
and non-universal ones for the third. Both these cases are well motivated from the 
GUT model building point of view. 

We explained how the remaining independent high-scale parameters are con- 
strained by the requirement of realistic electroweak symmetry breaking and a suf- 
ficiently large Higgs mass. We also briefiy compared with the related CMSSM and 
NUHM scenarios. Finally we presented a detailed parameter scan using Markov 
Chain Monte Carlo methods, highlighting the preferred ranges of parameters as 
well as correlations between them. 

Our analysis shows that models with degenerate Higgs mass matrix can be 
viable UV-completions of the MSSM for large ranges of gaugino and sfermion 
soft terms. They are, however, already strongly constrained by direct Higgs and 
SUSY searches, fiavour physics, and cosmology (as is the MSSM as a whole). 
In particular, the need to evade the LEP Higgs mass bound leads to preferred 
sparticle masses in the TeV range. This implies large finetuning in obtaining the 
correct electroweak scale. Most of the parameter points we found are fine-tuned 
on the sub-percent level, which of course refiects nothing but the well-known little 
hierarchy problem of the MSSM. Another stringent constraint arises from the dark 
matter relic density: In the models we considered, the neutralino relic density is 
generically larger than the observed value, so rather special parameter values are 
necessary in order to enhance the neutralino annihilation cross section. 

Nevertheless, we find Qh^ c^ 0.1 over a large part of the parameter space. 
This is mainly due to Higgs funnel annihilation, a large Xi higgsino fraction, or 
coannihilation with sleptons. In the case of universal sfermion soft terms, the 
Higgs funnel is clearly the most important process. Here it is worth noting that 
the shapes of the ID posterior probability distributions are more or less generated 
by just demanding correct EWSB, with the other constraints adding little to the 
shapes. In other words, the EWSB condition already selects the parameters such 
that most of the low energy observables are of roughly the correct magnitude, with 
exception of the relic density. It is then mainly the relic density constraint that 
helps shape the likelihood maps, and this reshaping can be understood in terms 
of the different (co-) annihilation channel contributions. The global features of the 
probability distributions are also quite robust against the fine-tuning prior. 
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Most of the parameter space lies within the reach of LHC at 14 TeV. In the 
region which is most difficult for the LHC to access, the LSP is higgsino-like 
and spin-independent direct dark matter detection experiments should soon see 
a signal. Should the MSSM with degenerate Higgs mass matrix be realized in 
nature, it will therefore almost certainly be observed within the next few years. 
This naturally raises the question of model discrimination: Can we look for a 
piece of experimental evidence pointing more or less uniquely to DHMM models? 
Unfortunately it seems to us that there is no such "smoking gun" signature for this 
kind of scenario. LHC may be able to exclude our models, but even if, conversely, 
an MSSM spectrum compatible with DHMM was found, it would need a future 
linear collider to accurately measure the sparticle masses and make a bottom-up 
reconstruction of the GUT-scale structure feasible. 
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